Coupling of Josephson current qubits using a connecting loop 
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We propose a coupling scheme for the three- Josephson junction qubits which uses a connecting 
loop, but not mutual inductance. Present scheme offers the advantages of a large and tunable 
level splitting in implementing the controlled-NOT (CNOT) operation. We calculate the switching 
probabilities of the coupled qubits in the CNOT operations and demonstrate that present CNOT 
gate can meet the criteria for the fault-tolerant quantum computing. We obtain the coupling strength 
as a function of the coupling energy of the Josephson junction and the length of the connecting loop 
which varies with selecting two qubits from the scalable design. 

PACS numbers: 85.25.Dq, 03.67.Lx, 74.50.+r 



I. INTRODUCTION 

The persistent current qubit using superconducting 
loop with three Josephson junctions has been proposed 
as a promising candidate for quantum computer 1 ' 2,3,4 
owing to the advantages of relatively long decoherence 
time and possible scalability. The current qubit has 
shown quantum superposition 3 and, recently, coherent 
time evolution^ between two quantum states. Therefore, 
the feasibility of the current qubit as a practical quan- 
tum computer has increased recently. In order to realize 
a practical quantum computer, qubits must be coupled 
with scalability and any selected pair of qubits should 
perform a universal logic gate. One of universal logic 
gates is composed of the controlled-NOT (CNOT) oper- 
ations and the single-qubit rotations^*^ 

In this study we suggest a coupling scheme for CNOT 
operation and scalable design using the three-junction 
current qubit as our unit qubit whose feasibility as a sin- 
gle qubit has already been testediSi 3 ^ Here we propose 
a design for coupling two qubits using a connecting loop. 
Two qubits are connected by a superconducting loop with 
Josephson junctions which carries a persistent current by 
effective piercing flux. The effective flux comes from the 
phase differences across the Josephson junctions common 
to the qubit loop and the connecting loop. We calculate 
coupling strength between two qubits as a function of 
the coupling energy of Josephson junction in connecting 
loop and the length of connecting loop which varies with 
selecting qubits from the scalable circuit. 

Recent experiments on two coupled charge qubits 
have demonstrated coherent oscillation! and CNOT gate 
operation. 8 For the coupled three- Josephson junction 
qubits 9 and the coupled current-biased Josephson junc- 
tion qubits 1 ^ spectroscopy measurements have been per- 
formed, but the CNOT operation has not yet been 
demonstrated. It seems that, for the coupled three- 
Josephson junction qubits, the coupling strength between 
qubits is too weak to implement CNOT operation. Here 
we aim to propose a coupling scheme which is expected 
to have strong advantages in CNOT operation, since 
the coupling strength can be large enough to suppress 
the unwanted oscillation of the control qubit. In order 



to obtain quantitative results we calculate the switch- 
ing probabilities of the coupled qubit system in CNOT 
operation. For an appropriate range of parameters, tar- 
get qubit oscillates and shifts to another state while the 
unwanted oscillation of control qubit can be suppressed, 
which can be implemented within the fault-tolerant crite- 
ria of quantum computingA& Moreover present coupling 
scheme offers the tunable coupling strength and the se- 
lective coupling from the scalable design, which is crucial 
for implementation of practical quantum computing but 
is not provided by the previous coupling schemes^ 
We also show that the coupling strength does not dimin- 
ish rapidly as length of the connecting loop increases. 



II. HAMILTONIAN OF THE COUPLED QUBITS 

A pair of coupled qubits is drawn in Fig. 1 (a). Two 
qubits are coupled through a connecting loop. Using this 
coupled qubits as a basic unit, a possible design to achieve 
a scalable computing is presented in Fig. 1 (b). We ex- 
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FIG. 1: (a) A pair of coupled qubits. (b) A part of the scalable 
quantum computing device. One unit is composed of a three- 
junction qubit, a dc-SQUID, and a connecting loop. 
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press the direction of the penetrating flux of the qubits 
oppositely in Fig. 1 (a), because two selected qubits are 
connected in a twisted way as shown in Fig. 1 (b). The 
two Josephson junctions in the connecting loop of Fig. 1 
(a) denote the dc-SQUID's that enable tunable coupling 
and switching function for selective coupling. The wiring 
scheme connecting any pair of qubits is topologically the 
same. The interqubit magnetic crosstalk between the 
qubit islands can be neglected, since the induced flux of 
the three- Josephson junction loop is very weak and the 
qubit islands are separated far enough. In Fig. 1 (a) 
the weak magnetic interaction between the qubit loop 
and the connecting loop may exist even though the area 
of the connecting loop is very small. However, by man- 
ufacturing the connecting loop in such a way that the 
normal directions of the qubit loop and the connecting 
loop are perpendicular to each other, we can prevent this 
unwanted interaction. 

The fluxoid quantization condition for a thin super- 
conducting loop^i, 



$t + {m c /q c ) f v c ■ dl = n$ Q 



(1) 



with q c = 2e, m c — 2m e , v c the average velocity of 
Cooper pairs, $o = h/2e the superconducting unit flux 
quantum and <i>t the total flux, can be represented as the 
periodic boundary condition for the wave function of the 
Cooper pairs including the phase evolution accumulated 
over the circumference of the loop with three Josephson 
junctions as follows ^ISii 3 . 



k n l = 2nn + 2nft 



+ 4>2 + 03), 



(2) 



where k n is the wave vector of the Cooper pairs, n an 
integer, I the circumference of the loop, and ft = $f/$o- 
Here fa is the phase difference across the i-th Josephson 
junction in a qubit. 

Since the total flux is the sum of external and induced 
one, it is written as 



f f _l LsI t u 



(3) 



where / = ^cxt/^o and 7 = L s /Lk- Here L s is the self 
inductance and Lk = rn c l/An c q 2 the kinetic inductance 
of the loop, where n c is the density of Cooper pairs, A 
the cross section of the ring, and m c — 2m e the mass of 
a Cooper pair. The dimensionless parameter 7 can also 
be represented by the ratio of the energy scale such as 
7 = E K /E L with E K = $l/2L K and E L = <S> 2 J2L S . 
Here use of the current expression / = —n c Aq c hk n /m c 
has been made. Then k n l of Eq. (J5J| is reexpressed as 

(1 +l)k R (L)l = 

2vrn fl ( L) + 27t/r(l) - (</)r(l)i + 4>r{l)2 + <f>n(L)3), ( 4 ) 

where R(L) denotes the right (left) qubit and &R(n is 
k n l of the right (left) qubit. 



Since the value of 7 for the three-Josephson junction 
qubit loopii is very large such that 7 « 200, Eq. Q can 
be written as ^k R[L) l w 2-kti r{l) + 2irf R ( L) - (<t>R(L)x + 
§ R (l)i + §r(l)z) an d the boundary condition in Eq. |2"jl 
becomes 



2t™ + 2irf t - (0i + 4>2 + <fo) w 0. 



(5) 



Using the relation in Eq. @ we can obtain the relation 
for the current in the qubit loop given by 



R (L) 



$0/ , , 4>R(L)l +(t>R{L)2 + <t>R(L)3\ ,. 

-j- [ n R (L) + }R{L) 7^z ) • ( 6 ) 
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From this expression we can get for the effective potential 

U eS (<t>), 



u. 



R(L). 



cir 



^fl(L)l, <t>R{L)2i <t>R{L)3) 
3 

= Y J E Ji {l-cos(i>R {L)i )+E^ 1 



(7) 



where Eji — Ej R i — Ejh is the Josephson coupling 
• '" = Ift, L j/2L s can be represented as 



energy and E„ 
E R(L) = 
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In order to obtain the effective potential of connecting 
loop, we use the periodic boundary condition over the 
circumference of connecting loop with self inductance L' s 
and obtain the relation for the current in connecting loop 
such that 



1 = -LL ir 
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(9) 



where r is an integer. The superscript prime in this work 
means the connecting loop. The effective potential of 
connecting loop with the Josephson coupling energy E'j 
can be represented by 

u' cS {(l>' L i(t> R ARiALi) 

= E' n + E'j(l - cos^) + E'j(l - cos^), (10) 



where 
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The total effective potential, U^g, of the coupled qubit 
system is given by 

t/^(0) = U^(cf> LU <f> L2 ^ L3 ) + U^{ct> Rl ,cj>R2, fl3 ) 
+ C/^ ff (0' L ,^,0 m ,0 L1 ), (12) 



where <j> 



t>Ll, 4'L2,4>L3,4>R1,4>R2, 4>R3, <j>' L > $r) 
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The local minima of the total effective potential can 
be obtained by minimizing the effective potential U^g(<pi) 
with respect to fa, i.e., dU^g/dfa = 0, which results the 
current relations at the Josephson junctions as follows: 



/' + Ir 
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(13) 
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where i=2,3. 

The energy of connecting loop depends on current state 
of each qubit, i.e., diamagnetic or paramagnetic. We ex- 
press the diamagnetic persistent current state as | J,) and 
the paramagnetic state as | f ) using spin language. In 
this scheme one of three junctions in each qubit is in- 
cluded in the connecting loop as can be seen in Fig. 1 
(a). The phases of each qubit, <Pli and 4>ri, determine 
the energy of coupled qubits. Sign of these phases is posi- 
tive (negative) for diamagnetic (paramagnetic) persistent 
current. The value of ((f>m — 0li) in the periodic bound- 
ary condition of Eq. © for the connecting loop plays 
the role of effective external flux in the connecting loop. 
When one qubit is in the paramagnetic current state and 
the other in the diamagnetic state, i.e., | ||) or | f|), 
the phases 4>r\ and (p^i have opposite sign and induce a 
large value of (cj)Ri — 4>li)- In this case, the phase differ- 
ence fa L and 4>' R induced by this effective flux gives the 
large Josephson coupling energy in the connecting loop. 
However, when both qubits are in the same current state, 
i.e., | ||) or | ||), the phases 4*m and 4*li have the same 
sign and are canceled in {4>ri — <Pli)- Then the energy 
of connecting loop become very small. 

The charging energies of Josephson junctions give the 
kinetic part of Hamiltonian. The number of excess 
Cooper pair charges on Josephson junction iVj = Qi/q c 
is conjugate to the phase difference fa and the commu- 
tation relation [fa , iVj] = i gives the quantum phase fluc- 
tuations of the junction. Using the Josephson relation, 
Q t = C($o/2tt)0 1 , charging energy, E c = E,:<9?/ 2C *> 
becomes 




where C^^i and C" L ^ are the capacitances of the 
Josephson junctions of the qubit loop and the connect- 
ing loop, respectively. If we introduce the canonical mo- 
mentum Pi = Nih = —ihd/dfa and the effective mass 
Mij = ($o/2Tr) 2 CiSij, the Hamiltonian can be given by 

H=±P T -M- 1 -P + U^ifa), (18) 



which describes the dynamics of the particle with effec- 
tive mass M in the effective potential [/^(cp)?- 

For a given value of the external flux the effective po- 
tential of a three- Josephson junction qubit has a double 
well structure corresponding to the diamagnetic state, 

|), and the paramagnetic state, | f) at the potential 
minima. The diamagnetic state and the paramagnetic 
state have the energy difference depending on the value of 
the external flux. The charging energy Ec causes quan- 
tum tunneling through the barrier of the double well po- 
tential resulting the quantum mechanically superposed 
states. The Hamiltonian describing a qubit then becomes 
-ffqubit = ho<7 z + t a x , where 2h is the energy difference 
between two states, | |) and | |), and to is the tunneling 
amplitude. 

For the coupled two qubits, on the other hand, we have 
four states, | ||), | ||), | ||) and | ||), at the effective 
potential minima for a given value of (/L,/i?)- Solving 
coupled equations (|13|) — <|16|) in conjunction with Eqs. © 
and (jnj numerically, we obtain the local minima of the 
effective potential of the coupled qubits. We set 7 = 
200, L s = llpH and £:,7 4 =200GHz, which are the design 
parameters of the qubit of Ref. 3. For the connecting 
loop we set 7' = 7 and L' s = 5L S . Obtaining the effective 
potential minima, [/^((/^min), we get the energy levels 
at each potential well, 

E s = -Tk^ + U^OK,^), (19) 

where s denotes the four states, | ||), | ||), | |f) and 
ft) j and 4> Sim - m is the value of <p s at the potential min- 
imum corresponding to the state s. The ground state 
energy, (l/2)?ia; s , can be obtained in harmonic oscilla- 
tor approximation^ Since the plasma frequency of the 
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FIG. 2: Ground state energy diagram for coupled current 
qubits. A,B,C, and D denote the resonance lines and the ex- 
ternal flux difference between resonance lines, Sf g , originates 
from the coupling energy. We consider E'j = Ej. 
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Josephson junction of the three-Josephson junction loop, 
u>i = Eji/ (Qq/2-k) 2 C , is larger than the characteristic 
frequencies, i.e. the level spacing and the amplitude of 
the tunneling between two states | f) and | |), by or- 
der of two£ and the plasma frequency of the Josephson 
junctions in the connecting loop can be the same mag- 
nitude as that of the qubit loop, we can consider only 
the ground state energy. By calculating the energy E s in 
Eq. for a certain range of (/l,/^), we can obtain 

four energy planes. In Fig. El we show the lowest energy 
planes and the states of the coupled qubits corresponding 
to the energy levels. The finite value of Sf g comes from 
the coupling energy which is higher for different qubit 
state than the same qubit state. Similar diagrams are 
shown for coupled charge qubits^S^ 

The Hamiltonian for the coupled qubits can be written 
in the two-qubit basis, | ||), | ||), | IT), and | TT),£ from 
the calculated values of E s in Eq. (|19|l as follows, 



H, 



coup 




(20) 



where th(H) is the tunneling amplitude across the barrier 
of the double well potential of the left (right) qubit and 
can be calculated by various methodsi2ii£ These basis 
states we use do not directly correspond to the individual 
states without coupling, but to the coupled qubit states 
characterized by the current states of the two qubits and 
the connecting loop. A straightforward calculation leads 
another form of the Hamiltonian H rnur ,, 



-Hcoup = h L a z L (g) I + h R I •Sktr- Jcr£ <g> aj 
+ t L a x L ®I + t R I®a x R 



(21) 
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FIG. 3: Energy diagram of the current states of coupled 
qubits, when /l = fa — f and E'j = 0.02Ej. The ener- 
gies of the coupled qubits without tunneling are represented 
by solid lines. Dashed lines denote the eigenvalues of H coup 
in Eq. (1201 with th = tfl=lGHz. The energy of antiparallel 
spin state of coupled qubits shifts upward by amount of 2J 
compared with that of the uncoupled state shown as dotted 
line. 



with relations that E±± = hL + h R — J, Ef± — —hL + h R + 
J, E^ = hi, — h R + J and E^ = —h^ — h R — J. Here / 
is the 2x2 identity matrix and the coupling constant J 
can be represented as 



(22) 



which shows that J is the energy difference between the 
parallel spin states and the antiparallel spin states as 
shown in Fig. [31 

Figure [31 shows the energy levels of Fig. [21 along the 
symmetric line Jl — fa, where the eigenstates of the 
total Hamiltonian, H coup , are represented as dashed lines. 
The two-qubit coherent oscillations can be performed at 
the co-resonance point, fi = f R = 0.5. The coupled 
qubits are taken into the co-resonance point from the 
initial position far away from it and stay for a finite time 
during which temporal evolution occurs^ 



III. CONTROLLED-NOT OPERATIONS 

We here consider the CNOT operations using evolu- 
tions at resonance lines 8 in Fig. [21 and show that the 
large value of J in present coupling scheme is able to give 
accurate CNOT operations. In Fig. 0]the energy levels 
around lines A and B are shown. Initially the coupled 
qubits are located far away from the resonance lines, say 
Ir = Jl — 0.457, being represented by the linear combi- 
nation of the states, | 11), \ T4): I IT) an d | TT)- When we 
adjust f R to the resonance line B in Fig. 01 the evolution 
from | Tl) to | TT) takes place, and vice versa. How- 
ever the states, | ||) and | IT), do not respond to this 
operation. This discriminating operation enables one to 
implement the CNOT gate using the left qubit as a con- 
trol qubit and the right qubit as a target qubit. The 
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FIG. 4: Energy diagram near the resonance lines A and B 
when /l = 0.48 and E'j — O.lEj. Solid lines represent the 
energies of the coupled qubits without tunneling and dashed 
lines are for tr, = tfl=lGHz. A and B show the values of fn 
where the point (/l,/a) is on the resonance lines A or B. 



5 




0.0 0.5 1.0 1.5 2.0 



nt/n 

FIG. 5: Time evolutions of the coupled qubits with tn = 
£ L =0.8GHz for (a) E'j = Ej,f R = 0.457 and f L = 0.4837 
and (b) E'j = 0.005Ej,f R = 0.4995 and f L = 0.4793. The 
solid (dashed) lines display the probabilities that the target 
(control) qubits occupy the corresponding states. 

resonance lines C and D in Fig. [21 are also available, if 
we want to use the right qubit as a control qubit. 

In order to demonstrate the CNOT gate we study the 
time evolutions of the coupled qubits. We start with the 
two qubit state written by 

=ci\ ||) + c 2 | U) + c 3 | IT) + c 4 | TT>. (23) 

Considering the CNOT gate operation at resonance line 
B in Fig. we set the initial condition, c\ = Ci = 1 
and C3 = C4 = 0, at t ~ 0. The time evolution of this 
state is obtained by the Schrodinger equation, H^>(t) — 
ihd^(t)/dt, with the Hamiltonian in Eq. J20|). In Fig. 
|5] (a) we show the time evolutions of \ci\ 2 for the case 



E'j = Ej and t L = t fl =0.004i;j=0.8GHz. At time fit w 
7r with fl = 2tji = 2<£ we can see that the state | fj.) 
evolves into the state | ft), while the states, | J, J,) and 
| jt) do not respond this operation. In order to invoke 
the spin flip between the states | J, J,) and | ||) remaining 
the states | f I) and | ff) unchanged, the resonance line A 
can be used for the similar process. Therefore the present 
scheme completes the CNOT gate operation. 

For this discriminating operation the resonance lines 
A and B should be kept far away from each other. When 
the target qubit evolves in the CNOT operation, the con- 
trol qubit must stay in the initial state. However, if the 
resonance lines are too close, the unwanted oscillation of 
the control qubit may take place. Thus the design pa- 
rameters of the coupling scheme should be set for the 
value of 6f g to be large enough to exhibit high perfor- 
mance of the CNOT gate. In Fig. 0] when /r is adjusted 
to line B, the amplitude of the unwanted oscillation be- 
tween two states, | J. J.) and | jf), can be suppressed by 
the large value of Sf g (or J)£ We obtain the values of 5f g 
and J as a function of E'j as shown in Fig. EJ We also 
consider the small Sf g case when E'j = 0.005i?j=lGHz 
resulting J=0.49GHz to investigate how the small cou- 
pling strength deteriorates the CNOT gate performance. 
In Fig. [3] (b) we show the time evolution of the this cou- 
pled qubits with the same initial condition as Fig. [5] (a) 
and tL = ti?=0.8GHz, where we can see that the con- 
trol qubit also evolves and does not stay in the initial 
state when the target qubit flips at fit « n. Thus we 
see that the CNOT gate cannot be accomplished in this 
parameter regime. 

Present coupling scheme can achieve large coupling 
strength by using the phases, not flux, of the phase qubit. 
With the phase differences, 4>li and 4>m, of the Joseph- 
son junctions in Fig. ^ ( a ), the effective flux has the 
value that f cS = (<j> R1 - cf) L1 )/2ir « 0.32 for f L « 0.5 




1.0000 



0.10 



CO 



0.05 



0.00 



FIG. 6: Plots of the coupling strength J and the external 
flux difference Sf g versus the Josephson coupling energy in 
the connecting loop. Inset shows the the coupling strength J 
versus interqubit length I' /I for E'j = 2Ej. Here the coupling 
strength J is calculated at the co-resonance point /l = fn — 
0.5. 
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FIG. 7: Switching probabilities of the state | jj.) to the state 
I |t) °f the coupled qubit system when E'j — Ej and /l = 
0.457 for various tunneling amplitudes tL = t R — to. Inset 
shows the probabilities that the state | J. J.) does not respond 
to this operation and remains unchanged. 
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and fn s=s 0.5. The value of coupling strength J can be 
about O(Ej) with i?j=200GHz by varying the value of 
E'j as shown in Fig. EJ while it is just about 0.5GHz in 
the experiment using the mutual inductance- since the 
induced flux of three- Josephson junction loop is quite 
smallii such that / in d = $ind/$0 ~ 0.002. 

The remaining problem is about the accuracy of 
CNOT operations. By deliberately choosing the pa- 
rameter regimes for the implementation of CNOT gate, 
we can satisfy the fault-tolerant criteria of quantum 
computation 5 which is about 10~ 4 . We calculate the 
switching probabilities for the case in Fig. [3] (a) when 
ilt ps 7r. In Fig. the switching probabilities for vari- 
ous tunneling amplitudes are shown as a function of the 
external fluxes, fu, with the fixed value of Jl = 0.457 
and E'j = Ej. We see in this figure that the fidelity is 
higher for lower tunneling amplitude. For higher tunnel- 
ing amplitude, the evolution at a resonance line can be 
affected by the presence of the other nearby resonance 
lines. Actually the external fluxes in Fig. [7| correspond 
to the point in the sector of resonance line B between 
the lines C and D where the evolution of the left qubit 
occurs (Fig. [2J. The large tunneling amplitude ti, may 
induce the unwanted oscillation of the controll qubit and 
reduce the fidelity for the state of the target qubit. The 
slight difference of pick points in Fig. [7|is also attributed 
to this aspect. We note that, although the probabilities 
are not so high as in the low to cases, the high to cases 
have the advantage that the operating ranges of fa are 
rather broad than the low to cases. In the inset of Fig. 
|7|we show the probabilities that the the state | J, J,) stays 
in the initial state. 

Many two qubit operations need the controllable cou- 
pling strength in implementation. However, the re- 
cent experiment on coupled charge qubits^ and on 
the coupled current qubits^ cannot control the cou- 
pling strength. For charge qubits a controllable cou- 
pling scheme using variable electric transformer has been 
proposed theoretically 14 Present coupling scheme allows 



the tunable coupling strength for superconducting cur- 
rent qubit. The Josephson junction in connecting loop is 
a dc-SQUID which is composed of two Josephson junc- 
tions with coupling energy Ej x . Total Josephson cou- 
pling energy E'j — 2Ej x cos(nf x ) is controllable by the 
finite flux f x piercing into the dc-SQUID. We show the 
coupling strength of the coupled qubits J in Eq. IL'l'l) as 
a function of E'j in Fig. 

In addition, we study the dependence of the coupling 
strength J on the length of the connecting loop, since 
the length of each connecting loop I' varies with selective 
coupling in the scalable design of Fig. Q (b). We plot 
J as a function of V jl in the inset of Fig. HO where J 
remains almost constant as the value of V jl increases. 
This behavior show that present scheme is appropriate 
for large scale integration which requires accurate two 
qubit operations even for very remote qubits. 

IV. SUMMARY 

We have proposed in this work a scheme of qubit cou- 
pling for the persistent current qubits using a connecting 
loop with Josephson junctions. The Hamiltonian of the 
coupled qubit system is obtained by calculating the en- 
ergy levels with the effective potential energy and the 
kinetic energy originated from the charging energy of the 
capacitance of the Josephson junctions in qubits and in- 
troducing the tunneling through the barrier of the double 
well potential. The time resolved dynamics shows that 
the accurate CNOT gate operations can be performed 
due to the large coupling strength which can be as large 
as O(Ej). Moreover the coupling strength is controllable 
and has a weak dependence on the length of connecting 
loop of two coupled qubits selected from the scalable de- 
sign. 
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